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I present some new results regarding confinement as it appears in Coulomb gauge. It is found that: 
i) a recently proposed Yang-Mills vacuum wavefunctional in temporal gauge and 2+1 dimensions 
yields a Coulomb-gauge ghost propagator and linear Coulomb potential in good agreement with 
lattice Monte Carlo results; ii) adding a few constituent gluons to heavy quark-antiquark states 
brings the interaction energy much closer to that of the static quark potential, and suggests the 
beginnings of gluon-chain formation at roughly one fermi; iii) a perturbative approach to Faddeev- 
Popov eigenvalues indicates that the zero eigenvalue at the Gribov horizon may occur either at, or 
away from, p = 0, depending on the gauge choice and spacetime dimension. This last result may 
be relevant to the qualitatively different infrared behavior of the ghost propagator in Coulomb and 
Landau gauges. 
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1. Introduction 

In this contribution I would like to report on three recent results concerning the non-perturbative 
behavior, in particular the confinement property, of Yang-Mills theory in Coulomb gauge. The first 
result, regarding the Yang-Mills vacuum wavefunctional and the Coulomb potential that can be 
derived from it, and the second, concerned with the energetics of physical states containing con- 
stituent gluons and a static quark-antiquark pair, were obtained in collaboration with Stefan Olejnfk. 
The latter work is presented in more detail in ref. [jj]]. In the last section I will outline another new 
result, concerning the spectra of the Faddeev-Popov operator at the Gribov horizon in less than four 
dimensions, in Coulomb and Landau gauges. 



2. The Vacuum Wavefunctional and the Coulomb Potential 

A long time ago it was argued [Q] that the pure Yang-Mills vacuum wavefunctional, in temporal 
gauge, has the following form at large distance scales: 1 



% ff \A] «exp 



d 3 x FfjFfi 



(2.1) 



This vacuum state has the property of dimensional reduction, in the sense that computation of 
a large spacelike loop in 3+1 dimensions reduces to the calculation of a Wilson loop in Yang- 

(2) 

Mills theory in 3 Euclidean dimensions. Assuming the vacuum state *J/q of the 2+1 dimensional 
theory also has the dimensional reduction property, the calculation of a planar spacelike loop in 
3+1 dimensions reduces to a calculation in D=2 dimensions, i.e. 



W(C) = (Tr[U(C)]) D = 4 
~ (Tr[U(C)]) D = 3 
~ (Tr[U(C)]) D = 2 



(^ 3) |Tr[f/(C)]|^ 3) ) 
(¥V\Tr[U(C)]\¥V) 



(2.2) 



In D=2 dimensions the Wilson loop can be calculated analytically, and we know there is an area-law 
falloff. Dimensional reduction therefore implies confinement. 

Support for the dimensional reduction form of the vacuum wavefunctional comes from strong 
coupling expansions, and from numerical simulations. At strong couplings, in temporal gauge, the 
Yang-Mills ground-state has the form [@] 



W [U]=exp 



-l^TrlUUUWj + Oig 6 

° plaq 



(2.3) 



which clearly has the dimensional reduction property. At weaker couplings, a numerical approach 
is available We consider a modified lattice Monte Carlo simulation in temporal gauge, in which 
the configurations on the t = time-slice are restricted to belong to a finite set {U^ m \x),m = 
1,2, ...,M}. After each sweep through the rest of the lattice, the configuration at t = is chosen 



A similar proposal was made by Halpern in 2+1 dimensions pfl 
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from among the given set via the Metropolis algorithm. Let N m be the number of times the ra-th 
configuration is selected. Then it is easy to show that in the limit that the number of sweeps Nj — > 00 



lim 



(2.4) 



For "large-scale" configurations, e.g. non-abelian constant lattices and long-wavelength plane 
waves, the numerical results agree perfectly with the dimensional reduction form ( |2~I| ) 

Of course, the true vacuum can't be simply the dimensional reduction form; that would give 
incorrect results at short distances and high frequencies. Stefan Olejnfk and I have proposed, in ref. 
[Hi, that the ground state solution in D=2+l dimensions, in temporal gauge, is approximated by 



^o[A]=exp 



I /WW 



1 



ab 



yJ-D 2 -h ) + m 2 



B»(y) 



(2.5) 



where B a = F" 2 , D 2 is the covariant Laplacian in adjoint representation, Xq is the lowest eigenvalue 
of — D 2 , and m is a free parameter chosen from a fit to the string tension. The motivations and 
successes of this wavefunctional, and the method which was developed to simulate it numerically, 
are described in the cited reference (and in Stefan Olejnfk's contribution to this meeting [^]). Now 
the ground state wavefunctional in Coulomb gauge is simply the restriction of the temporal gauge 
wavefunctional to configurations satisfying the Coulomb gauge condition, i.e. *Fq°"'[A] = ^[A] 
with V • A = 0. If we denote by g the gauge transformation taking a configuration A to Coulomb 
gauge, then it is not hard to see that for any observable Q [A], 



(% oul \Q[A}\% oul ) 



(Vo|fi[poA]|Vo 



(2.6) 



Then (Q) can be evaluated by generating lattices with probability weighting *Fq, transforming these 
lattices to Coulomb gauge, and evaluating the observable. Lattice configurations which are gener- 



ated stochastically from the (latticized version of) *Po in eq. ( |2.5| ) will be referred to as "recursion 
lattices". The result for (Q) derived from recursion lattices can be compared with that obtained by 
generating lattices (which we refer to as "MC lattices") by the usual Monte Carlo procedure with 
a Wilson action, transforming to Coulomb gauge, and again evaluating the observable. This latter 
procedure is the standard method for obtaining the expectation values of observables in Coulomb 
gauge. 

We have computed both the ghost propagator G(R) and the Coulomb potential Vc(R), defined 

by 



GR) 



1 



VD 



, V C (R) 



1 



VD 



(-V 2 



1 



VD 



(2.7) 



x 57 \ v 7 xy/ 

The result for the ghost propagator, obtained from both recursion and MC lattices, is shown in 
Fig. [T[ The Coulomb potential is very sensitive to "exceptional" configurations with very small 
eigenvalues of the Faddeev-Popov operator —V • D; these lead to huge errorbars. To compare 
recursion and MC results, we impose cuts on the data, throwing away these rare configurations. 
The result is shown in Fig. |2[ Obviously, the results for the ghost propagator and Coulomb potential 
obtained from our proposed wavefunctional closely agree with those obtained by standard methods. 
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Figure 1: The Coulomb ghost propagator evaluated on both recursion and MC lattices, at lattice coupling 
P = 6 in D = 3 dimensions. 
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Figure 2: Coulomb potential evaluated from data sets with a cut on lattices with V (0) < —20, for p = 6 and 
lattice extension L = 24. Results (with the same cut) are shown for both MC and recursion lattices. 



3. Constituent Gluons and the Gluon Chain Model 

The color Coulomb potential is known to be linear from computer simulations [§], but there 
are (at least) two serious difficulties in claiming that the Coulomb potential explains confinement. 
First of all, the Coulomb string tension a c is about three times larger than the asymptotic string 
tension a. Secondly, since the Coulomb force is essentially a one-gluon exchange effect, there are 
inevitably long-range dipole forces, which would result in long-range van der Waals forces among 
hadrons. This latter problem is generic to any model of confinement based on ladder diagram 
exchanges, and the difficulty can be traced to the absence of a flux tube. This raises the question: 
given that we already have a confining Coulomb potential, how and why does a flux tube form in 
Coulomb gauge? 

Let's begin with the fact that the Coulomb potential Vc(R) is the interaction energy of a certain 
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physical state containing two static charges, namely 

%«=r(o)<f(W) (3.i) 

But this state is not necessarily the minimum energy state with two static charges, so the question 
is whether (and by how much) we can bring down the interaction energy, by adding some "con- 
stituent" gluons corresponding to acting on the vacuum with gluon field operators. Schematically, 
these produce physical states of the form 

= r(0){c + cjA + c 2 AA + ...y h q b (R)y (3.2) 
According to a scenario known as the "gluon chain model" (cf. ref. [fl and references therein, and 



ref. [10]) the QCD color electric flux tube can be thought of, at least in certain gauges, as a chain 
of constituent gluons. The idea is that as a quark-antiquark pair separate, the electric field energy 
rises until a point is reached where it becomes energetically favorable to have a gluon in between 
the two quarks, thereby reducing the effective color charge separation. As the quarks continue to 
separate, this process is repeated, and eventually the state resembles a chain of constituent gluons 
(Fig. |3|), with the ordering of gluons in the matrix product of color indices correlated with the spatial 
ordering of gluons between the color charges. This scenario has a number of attractive features. 
A gluon chain has string-like properties (e.g. a Luscher term), Casimir scaling is obvious at large 
N, and the model is also consistent with N-ality dependence, in that string-breaking is natural for 
quark-antiquark sources in higher-dimensional color representations, at sufficient separations 

i I I • - 

q S luon gluon gluon q 



- 0----0 • • 

^ q gluon q 

time 

- o---m - 

q q 
Figure 3: The gluon chain model 

We would like to know if constituent gluons do indeed lower the interaction energy, for quark- 
antiquark separations that can be achieved in numerical lattice simulations. Let us define, on the 
lattice, the rescaled transfer matrix 

T = exp[-(H-E )a] (3.3) 

where Eq is the energy of the ground state, a is the lattice spacing, and exp(-Ha) is the usual 
transfer matrix. Ideally, we would like to diagonalize T in the subspace of states containing two 
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static charges. In practice, we must diagonalize in a finite M-dimensional subspace. Let 

\k)=r(x)Q a k b q b (y)\^o) k=\,2,-,M (3.4) 

where the Qu are functionals of the link variables. In general the {\k)} are not orthogonal. We then 
use standard lattice Monte Carlo methods to compute the following quantities 



O mn = (m\n) 

(3.5) 

t mn = (m\T\n) 

= {^Tv[Qi(t + l)ul(x,t)Q n (t)U Q (y,t)\) (3.6) 

From these quantities we can construct, via the Gram-Schmidt procedure, an orthonormal set of 
states {|cpi)}> and and also derive the matrix elements 7], = ((p,|r|(p ; -). The next step is to numeri- 
cally diagonalize the T-matrix in this finite basis, and then take 

V(R) = -\og(X max ) (3.7) 

where X mffiC is the largest eigenvalue of the r-matrix, as our estimate of the static quark potential. 

We build the <2/t operators out of "smoothed" A-fields. Define, in the usual way, A |U (x) = 
(U^(x) — U^(x))/2i. These fields are Fourier transformed, and the high-momentum components in 
directions transverse to direction e 7 (direction of the xj axis) are exponentially suppressed accord- 
ing to 

A,-(k,/) -»• exp -p(k 2 -kj) Ai(k,t) 

-> exp[-pki]A;(k,?) (3.8) 

where p is a variational parameter. Transforming back to position space, we denote the result- 
ing "transverse-smoothed" operator Ai(x,t,j) as the A-field smoothed in directions transverse to 
direction ey. We also define B,-(x,f) = 1 — ^Tr[Ui(x,t)] and smooth in the same way to obtain the 
transverse-smoothed operator Bi(x,t,j) . The Q operators are then defined in terms of the A,-(x,f, j) 
and Bj(x,t,j), for an antiquark at site xo and a quark at site xo + Rej 

Qi(t) = 1 2 

R-l 

02 (0 = J^Aj(x + nej,t,j) 

n=0 

R+l R+l 

23(0= £ L M x o+ ne j^j) A j( x o+ n ' e j^j) 

n=—2 n'=n 
R+2 R+2 

04 (0 = L L J^A i (xo + nej,t,j)A i (xo+n / e j ,t,j) 

n=—2 n'=n i^j 
R-\ 

Q 5 (t) = £b 7 (x + nej,t, 1)1 2 

n=0 

R-\ _ 

0e(O = ££Bi(xo + ne 7 -,?,y)l 2 (3.9) 
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where we have also defined, for i ^ j,A,(x,t,j) = A (A; (x, t , j) + A,- (x — e, , t , j) ) . From these opera- 
tors we construct and diagonalize the rescaled transfer matrix, as described above, in a truncated ba- 
sis of six states. We choose the variational parameter p which maximizes, at each quark-antiquark 
separation, the largest eigenvalue X max of T. Denote the corresponding eigenmode 

6 

\y(R)) = |(p„) (3.10) 

n=\ 

Then the fraction of the norm of \|/ due to zero, one, and two or more constituent gluons is 
|°i| > | «2 1 . and 1 — \a\ | 2 — \a2\ 2 respectively. 
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Figure 4: The color Coulomb potential V cou l {R), the "gluon-chain" potential V (R) derived from the minimal- 
energy variational state, and the "true" static quark potential V trU e(R) obtained by standard methods. Results 
are shown at lattice coupling p = 2.4. Continuous lines are from a fit of data points. 



gluon content in the chain state 



1 



' ;:■.>:. 


°" •-? V, X. , 






- Og, p=2.2 + 
Og, (3=2.3 — x— 
Og, |3=2.4 

" 1g, p=2.2 □ 
1g, p=2.3 

- 1g, p=2.4 - -e - 
2g, p=2.2 • 
2g, p=2.3 

- 2g, p=2.4 « 


• ; .'::«k 

S#S!S» 

jSS"" 


















4 .. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 



R (fm) 



Figure 5: Zero, one, and two-gluon content (fraction of the norm of the variational state) vs. quark separation 
R in fermis, at P = 2.2,2.3,2.4. 

Figure ^ shows our result for V(R) at (3 = 2.4, together with the corresponding Coulomb 
potential (from — log(7\i)), and the static quark potential computed by the usual methods. We 
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Figure 6: Sensitivity of the Coulomb potential V cou i(R) (solid symbols), and insensitivity of the chain poten- 
tial V c hain(R) (open symbols), to lattice volume. Data is for the gauge coupling p = 2.4, and lattice volumes 
L 4 = 12 4 , 16 4 ,22 4 . Quark-antiquark separation R is in lattice units. 



see that i) constituent gluon operators bring the potential down from the Coulomb value to a value 
much closer to the usual static quark potential, and ii) the addition of constituent gluons does not 
affect linearity; the potential V(R) of the minimal-energy variational state is still rising linearly. 
Fig. |5| shows the gluon content of the minimal energy state, where we see that at separations of 
about one fermi the variational state contains an equal admixture of zero and one-constituent gluon 
states. Thus, at one fermi separation, we may be seeing the beginnings of a gluon-chain structure. 
This figure shows results obtained at (3 = 2.2,2.3,2.4, and is therefore also a test of scaling. 

The dipole problem has already been mentioned. The color Coulomb field is not expected to be 
collimated into a flux tube, and this means that there should be strong sensitivity to lattice volume, 
on a lattice of spatial extension L, for quark-antiquark separations close to R = L/2. The reason 
is that for separations of that size, the finite volume cuts off a region where the field energy is still 
significant. On the other hand, if the field energy were collimated into a flux tube of diameter d, 
and if L > > d, then there would not be a similar sensitivity to the finite volume. Figure ^ displays 
the Coulombic and the variational ("chain") potentials, computed on lattice volumes 12 4 , 16 4 ,22 4 . 
The results for variational, constituent-gluon states seem to be insensitive to lattice size, in contrast 
to the Coulomb potential. This may be a hint that the dipole problem is much less severe for the 
constituent-gluon states. 



4. Faddeev-Popov spectrum at the Gribov horizon 

"Coulomb confinement" means that the Coulomb energy of an isolated color charge is infinite. 
An important question is whether this property follows from dominance, in the Coulomb gauge 
functional integral, of configurations lying on or near the first Gribov horizon [11]. The first Gribov 
horizon is characterized by the existence of one non-trivial zero mode of the Faddeev-Popov (F-P) 
operator, and it is thought, since the ghost propagator is the inverse of the F-P operator, that an 
enhanced density of near-zero modes would result in a ghost dressing function which is singular at 
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p = 0. It is a bit of a puzzle why this should be true in Coulomb gauge but not also in Landau gauge, 
where numerical simulations strongly suggest that the ghost dressing function is non-sing ular @]. 
The question I would like to address is whether the near-zero modes associated with the first Gribov 
horizon necessarily correspond to singular behavior at p = 0. 

Ghost and gluon propagators have been studied extensively via Dyson-Schwinger equations; 2 
here I would like to investigate the spectrum of the F-P operator via ordinary second-order per- 
turbation theory, with non-perturbative information entering via an ansatz for the transverse gluon 
propagator. I work in less than four spacetime dimensions, to avoid the complications of renormal- 
ization. The F-P operator is 

M ac = _5^v 2 -g/ flfcc, Af(x)3 ; - 

= K + gK 1 (4.1) 
with eigenvalue equation M ab ($ b p = XpCp^, and 

(l p ) = 4 0) + (A4 1 ') + (AXf) + ... (4.2) 

with 4°^ = P 2 at zeroth order, and (A)S P , ) = at first order in g. The transverse gluon propagator 
is required for the second-order result, and I use Dfj(q) = 5 ab (dij — qiqj/q 2 )D{q) with the ansatz 
for the dressing function 



D(q) 



1 Coulomb gauge 



2y/q 2 +m 2 + a /q a 

(4.3) 



1 Landau gauge 



q 2 +m 2 + a /q a 

The second-order result in d + 1 spacetime dimensions (Coulomb gauge), or d spacetime dimen- 
sions (Landau gauge), is 

(X p )=p 2 (l-g 2 R d I\p,m,a}) (4-4) 
where Rj is a dimension-dependent constant of 0(1), and 

rn/2 

I\p,m,a} = / d& sin^ 2 6(l-cos 2 e) 
Jo 

x| f dc i 7;\D(4pcosQ + q)+D(q)] 

\Jo 4 + 2/?cos6 L v F H> WJ 

+ / d 4 o r[D(4/>cos6-<7)-D(<7)] (4.5) 

Jo q — 2pcos9 J 



with D{q) = q d ~ 2 D{q). Eq. (4.4) is somewhat reminiscent of the Dyson-Schwinger equation for the 
ghost propagator. However, this is an equation for the F-P eigenvalue spectrum, not the propagator, 
and of course there is no claim that the equation is exact. In general, as p — > 0, 

I[p,m,a] = a[m,a] — b[m,a]p s + ... 

(Ip) = (1 - a[m,a])p 2 +b[m,a]p 2+s + ... (4.6) 



2 This approach has a large literature. There exist "scaling" solutions to the Dyson-Schwinger equation, with an in- 
frared singular behavior for the ghost dressing function (cf. the review by Fischer [|lj]), and there also exist "decoupling" 
solutions JI4] ] in which the ghost dressing function is finite. At present, the numerical data in Landau gauge appears to 
favor the decoupling solution 112]. 
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Then in Coulomb gauge, for fixed a, I find three cases: 

1) m < m c , a[m,a] > 1, (X p ) ~ — p 2 



2) m > m c , a in.u. < 1, (X p ) ~ +p 2 



(4.7) 



3) m = m c , a[m,a] = 1, (X p ) ~ +-/? 2+ - ? 

When m is below some critical mass parameter m c , there is an interval of negative eigenvalues 
which begins at p = 0. An example at a = 1, Coulomb gauge in 2+1 dimensions, is shown in Fig. 
7(a) , with dimensional quantities in units of g 2 . This case corresponds to the dominant gluon field 
configurations lying outside the Gribov region. As m increases, the interval of negative eigenvalues 
shrinks in size, until at some critical m = m c the interval shrinks to one point. This would corre- 
spond to the dominant field configurations lying right on the first Gribov horizon. At m > m c we 
have X p ~ p 2 , and the relevant gluon field configurations lie inside the Gribov region. The spectrum 



above, at, and below m c , at near-zero values p, is shown in Fig. 7(b), again for a = 1 and Coulomb 
gauge in 2+1 dimensions. 
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Figure 7: F-P spectra at a = 1. (a) m — 0.20 < m c . There is an interval of negative eigenvalues in the region 
< p < 0.009. (b) X p at low p, for m below, above, and equal to the critical value m c — 0.2228. 

The m = m c case is the one of interest, because this is where the near-zero F-P eigenvalues rise 
with a non-standard power p 2+s , leading to an enhanced eigenvalue density and (it can be shown) 
Coulomb confinement in 2+1 dimensions for any s > 0. For any fixed a we can calculate s and 
the critical gluon mass parameter m c . For example, at a = 1 it is found that m c = 0.223, and the 
exponent 2 + s is determined from a straight-line fit to a log-log plot of X p vs. p, as seen in Fig. ||. 
From this plot, we find X p ~ p 2 53 at small p and m = m c , going back to X p ~ p 2 at larger p. 

Landau gauge is a little different. In two spacetime dimensions the F-P Landau gauge spectrum 
is qualitatively similar to the Coulomb gauge result in three spacetime dimensions. However, in 
three spacetime dimensions, for any a > 0, it turns out that the interval of negative eigenvalues at 
m < m c is located away from p = 0, and at the critical m = m c the non-trivial zero eigenvalue is at 
some p > 0. This situation is shown in Fig. ^. 
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Figure 8: Log-log plot of the spectrum of the Faddeev-Popov operator (2+1 dimensional Coulomb gauge), 
for a = 1 at the critical m c = 0.223. A best fit at p < 1 yields k p = 1.21/j 2 ' 53 . 
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Figure 9: The low-lying F-P spectrum in Landau gauge, in D = 3 dimensions and a = 1, for gluon mass 
parameter m above, below, and equal to the critical value m c . 

While the spectrum of F-P eigenvalues does not translate directly into a prediction for the 
behavior of the ghost propagator (because the momentum behavior of the F-P eigenmodes must 
also be taken into account), what we do see is that two different scenarios for the behavior of the 
eigenmode spectrum at the Gribov horizon are possible. In the first scenario, the non-trivial zero 
eigenvalue is at p = 0, while in the other, it is located away from p = 0. It is natural to conjecture 
that the first scenario is associated with the confining properties of the ghost propagator in Coulomb 
gauge, while the second scenario has something to do with non-singular infrared behavior of the 
ghost dressing function in Landau gauge. 
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